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Abstract 

Charged and neutral, pion and kaon electric polarizabilities are extracted from lattice QCD using 
an ensemble of anisotropic gauge configurations with dynamical clover fermions. We utilize classical 
background fields to access the polarizabilities from two-point correlation functions. Uniform 
background fields are achieved by quantizing the electric field strength with the proper treatment 
of boundary flux. These external fields, however, are implemented only in the valence quark sector. 
A novel method to extract charge particle polarizabilities is successfully demonstrated for the first 
time. 
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I. INTRODUCTION 



A staple component of electrodynamics courses is the electric polarizability. Neutral 
materials immersed in electric fields polarize. At the atomic scale, electron clouds distort 
creating microscopic dipole moments that oriente opposite the applied field to minimize 
the energy. This simple principle accounts for dielectric properties of materials, a range 
of intermolecular forces, and properties of atoms and nuclei in applied fields. At the fem- 
toscale, hadrons too polarize in applied fields, but only against the strong chromodynamic 
interactions confining their electrically charged quarks into hadrons. 

Understanding properties of hadrons quantitatively is formidable. Quark and gluon in- 
teractions must be treated non-perturbatively for which lattice QCD has been developed, 
see |l[ f° r a review. Low-energy properties of hadrons, however, can be described using an 
effective theory of QCD, based upon treating pseudoscalar mesons as the Goldstone modes 
arising from spontaneous chiral symmetry breaking. A picture of hadrons emerges from 
chiral dynamics: that of a hadronic core surrounded by a pseudoscalar meson cloud. As 
some pseudoscalar mesons are charged, polarizabilities of hadrons encode the stiffness of 
the charged meson cloud (as well as that of the core). The form of pseudoscalar meson 
polarizabilities is consequently strongly constrained by chiral dynamics 0,0,0). Beyond the 
leading order, however, results depend on essentially unknown low-energy constants, which 
currently must be estimated in a model-dependent fashion. For the case of the charged pion, 
confrontation of these results with experiment has proven difficult, e.g. from the original 
measurement Q, to the most recent p], extracted results disagree with predictions made 
using chiral dynamics. New results with higher statistics and the first kaon results are 
anticipated from COMPASS at CERN Q. 

Lattice gauge theory simulations provide a first principles approach to determine hadronic 
polarizabilities from QCD, crucially test predictions from chiral dynamics, and confront 
experiment. Indeed the unknown low-energy constants of chiral perturbation theory can be 
determined by matching to lattice QCD computations. Furthermore, the ability to vary the 
quark mass allows one to directly explore the chiral behavior of observables, investigate the 
convergence properties of the perturbative expansion, and thereby test the predictions of 
the effective theory. The highly constrained form for hadronic polarizabilities within chiral 
perturbation theory leads to a stringent test of low-energy QCD dynamics. 

Electric polarizabilities of neutral hadrons have been calculated with lattice QCD using 
the quenched approximation at pion masses greater than 500 MeV jsl, 0]. There has also 
been a fully dynamical calculation of the neutron electric polarizability at a pion mass 
of 760 MeV [lOj. These calculations do not employ constant electric fields but attempt to 
mitigate effects from field gradients by imposing Dirichlet boundary conditions on the quark 
fields in the time and/or space directions. Such an approach leads to uncertainties that 
are difficult to quantify. In this work, we report on calculations of pseudoscalar meson 
polarizabilities using lattice QCD with dynamical configurations. A salient feature of our 
computation is that it utilizes a periodic lattice action with everywhere constant electric 
fields. Our calculations of meson polarizabilities are the first such to include effects from 
dynamical quarks. At this stage, however, we are restricted to electrically neutral sea quarks. 
Correcting for this malady would require at least an order of magnitude greater computing 
power. 1 Furthermore, we demonstrate for the first time how to extract charge particle 



Pseudoscalar meson polarizabilities first depend on sea quark charges at next-to-next-to-leading order 
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polarizabilities from lattice two-point correlation functions. 

We begin in Section [III by describing the implementation of constant external fields on 
a lattice. The Appendix considers the effect of non-uniform fields. Next in Section 
we detail how lattice two-point correlation functions can be utilized to extract the electric 
polarizabilities of both charged and neutral particles. Details of our lattice study are then 
presented in Section [TV] and summarized in a conclusion, Section [V] 



II. CONSTANT FIELDS ON A LATTICE 

To produce a constant electric field, £ = £z, we use the Euclidean space vector potential, 

A^x) = (0,0,-^X4,0), (1) 

where £ is a real-valued parameter. The analytic continuation £ — > —i£m produces a real- 
valued electric field in Minkowski space. Generally this continuation cannot be performed 
using numerical data because of non-perturbative effects, e.g. the Schwinger pair-creation 
mechanism [l3j is absent in Euclidean space. We are interested, however, solely in quantities 
that are perturbative in the external field strength, for which the naive continuation produces 



the correct Minkowski space physics, see [14j for explicit details 



To implement the background field on the lattice, we modify the SU (3) color gauge links, 
Ufj,(x), for each quark flavor by multiplying by the color-singlet Abelian links, Ufj, (2), for 
the external field, namely 

U,(x)^U,(x)Uf\x), (2) 



where lljf\x) = exp[iQAn(x)\, where Q is the quark electric charge, and A^(x) is given in 
Eq. ([T]). As this multiplication is carried out on pre-existing gauge configurations, the sea 
quarks remain electrically neutral. This approximation is imposed because of computational 
restrictions which will not be rectified in the near future without a significant increase in 
resources. 

The inclusion of the field via Eq. ([2]) does not lead to a constant electric field. On a torus, 



constant gauge fields require quantization [15|, [la, [17|. The basic argument is as follows 



With periodic boundary conditions, 2 the action is defined on a torus, which is a closed 
surface. For the field we wish to implement, the only plane with non- vanishing flux is the 
23-24 plane. The total area of the 23-24 plane is /3L, where L is the length of the 23-direction, 
and (3 is the length of the 24-direction. Because the torus is a closed surface, however, there 
can be no net flux through the 23-24 plane (modulo 27r), i.e. $ = Q£f3L = 2nn, with n as 
an integer. This leads to the 't Hooft quantization condition 

£ = ?EL. (3) 



in the chiral expansion . It is thus possible to extract physical information from simulations with 

vanishing sea quark charges by utilizing chiral perturbation theory. As the current study is restricted to 
one volume and one pion mass, we leave this investigation to future work. 
2 The argument applies equally well to the case of twisted boundary conditions on the matter fields ip{ x ) of 
the form: ip( x + Lxj) = e lSj ip(x), and analogously for the time direction. Dirichlet boundary conditions, 
on the other hand, inevitably lead to problems. 
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Here we have used the down quark electric charge, qd = —1/3 e, and note that the up quark 
will necessarily encounter properly quantized fields when Eq. is met because q u = —2qd- 
The argument presented for constant gauge fields applies to a continuous torus, and 



must be modified for a discrete torus, see e.g. [18|, ll9|, [20j. On a discrete torus, each of 
the elementary plaquettes must be identical with value: exp(iQ£) to arrive at the constant 
electric field £. With Eq. the plaquettes are identical in the bulk of the lattice but not 
at the boundary, where there are L plaquettes with differing flux. Each of these plaquettes 
wraps around from 24 = (3 — 1 to X4 = 0, with the common value: exp[iQ£(l — (3)]. This 
unwanted flux can be eliminated on L — 1 of the plaquettes by including additional transverse 
links, Ufi (x), at the boundary, 

U,{x) — > Up(x)Ujp{x)UjP±(x), (4) 

with Ufi (x) = expliQSfixsSfiiSx^p-i]. Now the field through every plaquette is 8, with 
only one exception: the plaquette at the far corner of the lattice, (£3,24) = [L — — 1), 
that wraps around to x 3 = 0, and £4 = 0. The value of this plaquette is: exp[iQ£(l — (3L)], 
which is identical to the plaquette in the bulk of the lattice provided 't Hooft's quantization 
condition, Eq. ([3]), is met. We have previously demonstrated the effects of using non- 
quantized field values, finding non- negligible shifts in particle spectra [2l|. We summarize 
our findings in the Appendix. In this work, we implement the external field using Eq. (j3j), 
and quantized values for the field strength £ in Eq. This choice corresponds to a 

completely periodic lattice gauge action, and thus corresponds to a field theory at a finite 
(but low) temperature in the continuum limit. 



III. CORRELATION FUNCTIONS 



For a neutral particle, it is straightforward to calculate the electric polarizability using 
standard lattice spectroscopy [s|. 3 One merely matches the long-time behavior of Euclidean 
two-point functions, g(x4, £), computed in QCD to the expectations of the effective hadronic 
theory to deduce the particle's energy. The lattice two-point function has the form 

g(x 4 ,£) = J2(m^mO)e, (5) 

X 

where is an interpolating field for the particle of interest (e.g. = d^s for the K°), and the 
subscript £ denotes that the correlation function is determined in the background electric 
field. This correlation function is matched onto the correlator G(x±,£), in the hadronic 
theory, 

G(x 4 , £) = Z(£)e~ E{£)xi + Z\£)e- E '^ )xi + . . . , (6) 

where the ellipsis represents exponentially suppressed contributions beyond the first excited 
state. The ground-state particle's energy, E(£), has a series expansion in the external field 



3 Strictly speaking this is only true at infinite volume. At finite volume, there are additional effects stemming 
from boundary conditions and the compact nature of the external gauge field 22, 2^|. As we employ only 
one lattice volume to demonstrate our methods, we neglect these additional corrections at this stage. 
Further analysis with multiple volumes and multiple pion masses is needed to control these systematics. 
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strength 



E{£) = M + Uita E £ 2 - j^fa EEE £ A + ... 



(7) 



where M is the particle's mass, a E its electric polarizability, and a E EE is a multiple electric 
dipole interaction strength. Here the ellipsis represents terms at higher order in the strength 
of the field. The sign of the polarizability term (quadratic Stark shift) is positive due to 
our treatment in Euclidean space. The amplitudes, Z(£) and Z'(S), also have expansions 
in even powers of £. As explained below, we are forced in our particular computations to 
consider contributions from excited states, shown in Eq. (J6]). The energy, E', of the first 
excited state has an analogous weak field expansion in terms of the mass M', polarizability 



etc. 



When charged particles are subjected to constant electric fields, we again match the 
lattice correlation function, g(x±,£), to the correlator calculated in the hadronic theory. 
With sufficiently weak fields, quarks and gluons will still hadronize into a tower of states of 
the same quantum numbers, specifically of the same charge. For times, £4, long compared 
to that set by the excited state mass, 24 ^> r' ~ 1/M', the excited state contributions to the 
two-point function will still be exponentially suppressed (albeit not a simple exponential). 
For times beyond r', we can assume the two-point correlation function will be dominated 
by the ground state. As this state is charged, the behavior of the correlation function will 
have a more complicated form than a simple exponential falloff with time. 

For a relativistic scalar particle of charge Q, consider the single-particle effective action 
in the hadronic theory. As the particle is composite, there are both Born and non-Born 
terms in the action. The non-Born terms account for non-minimal couplings of the field 
to the particle, such as polarizabilities. These couplings can be summed, as in the case of 
a neutral particle, into the energy, E(£), defined above. The Born couplings additionally 
must be summed to arrive at the charged particle two-point function. For the field specified 
by Eq. (pQ), the equations of motion for a scalar particle are exactly solvable, and lead to the 
two-point function in the hadronic theory 



G(xi,£) = Z{£)D(x A ,E{£),£) + Z'(£) D(x A ,E J (£),£) + ... 



(8) 



with the ellipsis representing contributions beyond the first excited state, and the relativistic 
propagator of a charged scalar, D{x^E{£),£\ given by [lij 



D(x 4 ,E{£),£) 



dsA 



Q£ 



27rsinh(<5^s) 



exp 



coth(Qts) 



(9) 



where E(£) can no longer be interpreted as the energy but remains given by Eq. (17j) . Classi- 
cally E(£) is the rest energy of the charged particle. The analogy to a classically accelerating 
particle occurs at the level of the particle's action. Interpreting the Euclidean time behavior 
of the particle's motion on a compact space in terms of acceleration proves difficult. For 
Q = 0, this propagator properly reduces to Eq. (jSJ). For sufficiently weak fields, or equiv- 
alently short times, the 0(£ 2 ) termjn the series expansion of the correlator reproduces 



the non- relativistic result derived in [24|]. Due to our particular anisotropic lattices, we in- 
clude contributions to the two-point function from the first excited state thereby stabilizing 
the extraction of ground state parameters. The quantum numbers of the excited state are 
identical to the ground state, i.e. Q' = Q. 
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IV. LATTICE RESULTS 



To demonstrate our method for extracting meson polarizabilities from lattice two-point 
functions, we have employed an ensemble of anisotropic gauge configurations with 2 + 1- 
flavors of dynamical clover fermions 25, 26[. The ensemble we use consists of 200 lattices of 



size 20 3 x 128. After an initial 1000 thermalization trajectories, the lattices were chosen from 
an ensemble of 7000 spaced either by 20 or 40 to minimize autocorrelations. The spatial 
lattice spacing of these configurations is a s = 0.123 fm 125, 261, with a non-perturbatively 



tuned anisotropy parameter of £ = a s /a t = 3.5, where a t is the temporal lattice spacing. The 
finer temporal spacing affords us the ability to better fit non-standard behavior for two-point 
correlation functions, and is critical for this analysis. On the ensemble, the renormalized 
strange mass is near the physical value, while the renormalized light quark mass leads to a 
pion mass of m n ~ 390 MeV. 

On each configuration, we compute at least 10 up quark propagators, 10 down quark 
propagators, and 10 strange quark propagators with random spatial source locations. Mul- 
tiple inversions were made efficient using the EigCG inverter {27|. Interpolating fields at 

on a 



the source are gen erated from gauge-covariantly Gaussian-smeared quark fields [28|, 129 



stout-smeared |30J gauge field in order to optimize the overlap onto the ground state. Inter- 
polating fields at the sink are constructed from local quark fields. Each propagator is located 
with source time at r src = 0. Randomization of the source time location, while improving 
the statistical sampling, would complicate the extraction of both charged and neutral me- 
son correlation functions, as two-point functions are no longer time-translationally invariant. 
For charged particles, the correlator in Eq. (jSJ) is explicitly a function of the sink time-slice 
and not simply a function of the source-sink separation (the full dependence on source time 
is given in [14]]), while for neutral particles the violation of time-translation invariance arises 
from volume effects. 

The external field was implemented using Eq. (J4J), and propagators were computed for 
nine values of the field strength, n, corresponding to the integer appearing in the quantization 
condition, Eq. Q. We use n = 0, which corresponds to a vanishing external field, as well 
as n = ±1, . . ., ±4. On our lattices, the expansion parameter governing the deformation of 
a hadron's pion cloud is given by [14j 

^4V = 0.18n 2 . (10) 
ml) 

From the size of this parameter, we anticipate the need to include terms beyond quadratic 
order in the electric field expansion of hadron energies. In our analysis, we include terms up 
to quartic order. Larger lattices will be required for better control over systematics relating 
to the electric field expansion of observables. 

Meson two-point functions were obtained for each source location on a given configuration. 
Individual results for the multiple source locations on each configuration were then source 
averaged. This procedure was carried out for each value of the external field. To satisfy 
invariance under parity transformations, under which £ — > —£, we took the geometric mean 
of correlators calculated at n, and — n on each configuration. This reduced the set of fields 
to five magnitudes corresponding to the integers, n = 0, . . . , 4. The ensemble of correlation 
functions was then used to generate 200 bootstrap ensembles. Fits to the bootstrapped 
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ensemble were performed as described below. 4 



A. Neutral pion and kaon 



Bootstrapped correlators for the neutral pion and kaon were obtained using the procedure 
described above. As we use standard spectroscopy to determine the polarizabilities for 
neutral particles, we handle these mesons first. To facilitate the discussion, we consider the 
standard effective mass, given by 

where g(t, 8) is the bootstrap ensemble-averaged correlator. Error bars on the effective mass 
are determined using the bootstrap ensemble. Effective mass plots for the neutral pion and 
neutral kaon are shown in Fig. [TJ For the neutral pion, so far we have only calculated the 
connected part of the correlation function. The effective mass plot should exhibit a plateau 
over a range of time when the ground state saturates the correlation function. The temporal 
extent of our anisotropic lattices is f3 — 128 a t = 4.5 fm, which is considerably smaller than 
typical isotropic lattices, where (3 > 7.5 fm. As one must wait long enough for the excited 
states to drop out, the pion and kaon effective masses never plateau (or barely exhibit a 
plateau) because of the backward propagating image from the time boundary. 

To extract the ground state properties, we fit the correlation function g(t,£) using the 
two-state form of G(t, 8) in Eq. ([6]) augmented to include a backwards propagating ground 
state, and backwards propagating excited state (i.e. a sum of two hyperbolic cosines). We 
use a correlated chi-squared analysis to fit the time-dependence of the bootstrap ensemble 
of correlators. To determine the fit window, we use black box methods comparing single 



and double effective masses, see [31l . I32L l33| for details on the latter. We found the same fit 
window, 8 < t < 29, could be used for a given particle for every value of the field strength. 
Alternate fits on the same window without backwards propagating states result in a 1.5% 
shift of the neutral pion energies, and a negligible shift of the kaon energies. 

As the parameters Z(S) and Z'{8) enter the fit function G(t, S) linearly, we utilize variable 
projection (see jUJ for references) to reduce the number of fit parameters from four down 
to two, namely just the energies E(£) and E'(S). 5 We perform these two-state fits on the 
entire bootstrap ensemble arriving at an ensemble of energies for each magnitude of the 
electric field £, in particular { t E i {£)} for the ground state, where i indexes the bootstrap 
sample, i — 1, . . . , N. As the ensembles of configurations for different field strengths are 
generated from the same underlying lattice configurations, correlations between the energies 
for different field strengths will be significant and it is important to account for these. On the 
bootstrap ensemble of energies, we perform electric-field correlated fits to the energy function 



4 We have performed multiple differing procedures to analyze the data, of which only one is described 
in detail in the text. Throughout we will comment on the alternate procedures. Most notably, fits have 
been performed using a jackknife procedure to determine uncertainties, and with a separate analysis of the 
positive and negative values of the field strength. Effects of correlations in the data have been investigated 
by blocking neighboring configurations, and consistent results have been obtained. 

5 We also analyze correlation functions by fitting the effective masses with two states. These three param- 
eter fits give consistent results. 
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FIG. 1: Effective mass plots for the (connected) neutral pion and neutral kaon. Values for n 
correspond to the magnitude of the quantized electric field £ , Eq. ([3|). Curved bands show fits to 
the correlation functions. For each plot, the band spans the fit window, and the width is set by 
the uncertainty in the ground state energy. The flat bands shown are ground state energies with 
uncertainty for each field value. On a lattice of infinite time extent, the effective masses should 
asymptote to these bands over long times. 
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TABLE I: Summary of fit results for neutral meson two-point functions for 8 < t < 29. Here tt 
refers to the connected part of the correlation function, the fine-structure constant is af s . = e 2 /4ir, 
and 1 — P is the integrated chi-squared. All quoted values are averages over the bootstrap ensemble, 
and are given in dimensionless lattice units. For the electric polarizability, a^ tt = aE^af.s.itdl) -1 
and the higher-order coupling, = oleee{^- <xj s a t a t ) _1 10~ 3 . The first half of the table 

summarizes the time-correlated fits to the energies in each field, while the second half summarizes 
the field-correlated fits. The fits I and II are described in the text. The second uncertainty on the 
polarizabilities is an estimate of the systematic due to the choice of fit window as explained in the 
text. 



E{£) given in Eq. (|7j). With the ensemble average energies denoted by £ (£) 
we minimize the correlated chi-squared, namely 



X 



Y J [E{S)-E{£)]C-^ I [T,{£')-E{£' 

££' 



(12) 



with the field-strength correlation matrix, Cs s', given by 



N 



a 



£,£' 



N — 



T J2 [E(£)-^(£)\ [£(£') -EiiS' 



(13) 



Because all three fit parameters, M, oe, and cIeee, enter the fit function E(£) linearly, the 
chi-squared minimization can be done analytically. Fits to the energy function are carried 
out on the bootstrap ensemble, resulting fit parameters are averaged, and the uncertainties 
from fitting and bootstrapping are added in quadrature. We perform two different field- 
correlated fits as follows: (I) a fit to all five field strengths using Eq. ((7|), (II) the same fit 
function but with the largest field strength excluded. Finally, to estimate the systematics due 
to the choice of fit window, we performed uncorrelated fits to the electric field dependence of 
meson energies determined on adjacent fit windows. We chose the nine fit windows obtained 
by varying the start and end times by one unit in either direction. On each time window, 
we determined the electric polarizability. The systematic uncertainty on due to the fit 
window is estimated as the standard deviation of the extracted over the various adjacent 
windows. Fit details and extracted parameters are tabulated in Table [H 
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FIG. 2: Plots of correlated fits to the electric field dependence of neutral meson energies. For each 
field strength, the bootstrap averaged energies are plotted with error bars reflecting the uncertainty 
from statistics and fitting. Fits I and II to the ^-dependence are also shown with the plotted bands 
reflecting the uncertainty in the parameters appearing in Eq. ([7]). 



From the extracted polarizabilities, we can investigate the electric field dependence of 
meson energies. This is done in Fig. [2] for the neutral pion and neutral kaon. For the con- 
nected part of the neutral pion, we see downward curvature of the energy with respect to 
increasing S, while for the neutral kaon the energy is comparatively quite flat. In physical 
units, the polarizabilities oF E , and , are not consistent with naive expectations. To at- 
tempt a qualitative explanation for the size of the ground state polarizabilities, we compare 
our results with predictions from chiral perturbation theory. The neutral pion electric polar- 



izability at one- loop is negative [3J, [35|. While this is surprising, the one- loop polarizability 
arises solely from the disconnected contraction between quark basis r\ u and % mesons [11] . 
Hence the negative sign owes to group theory weight of r\ u versus rja in the pion interpolating 
field, 7T ~ -^(jlu — %)■ As we have only calculated the connected part of the correlator, 

chiral perturbation theory suggests that d 1 ^ is an order of magnitude smaller than the naive 
expectation. While our result is of this magnitude, it is of the wrong sign (the average of 
r] u and rjd polarizabilities should be positive). This negative value could arise from volume 
effects, which are known to be non- vanishing at next-to-leading order in chiral perturbation 
theory j23|. For the neutral kaon polarizability, the one-loop chiral computation vanishes, 
even with electrically neutral sea quarks [l^ |. Our extracted neutral kaon polarizability, 
however, is smaller than typical two-loop contributions. Because the dominant volume cor- 
rections arise from pion loops, we expect the neutral pion and kaon volume effects to be 
of the same size. If the negative result for the connected 7r° is due to volume corrections, 
then the near vanishing result for the K° could be due to a near cancelation between the 
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n a t E(£) l-P 


n a t E{£) l-P 


7T U U.UOyi^J U.DD 

1 0.0702(6) 0.46 

2 0.0718(8) 0.61 

3 0.0733(16) 0.93 

4 0.0497(129) 0.97 


J\ U U.Ui)0i>[O) U.iU 

1 0.0979(4) 0.77 

2 0.0982(7) 0.78 

3 0.0958(10) 0.98 

4 0.0927(23) 0.97 


7T+ a t M a l % tt u X eee 1 - P 


K+ a t M a 1 !" a l ^ E l-P 


I 0.0692(2) 18(4)(6) 24(10) 0.30 

II 0.0692(2) 16(3)(3) 17(10) 0.64 


I 0.0971(2) 8(3)(1) 17(5) 0.03 

II 0.0969(2) 16(4)(3) 40(9) 0.23 



TABLE II: Summary of fit results for charged meson two-point functions for 8 < t < 29. Entries 
are as in Tabled 

polarizability and the volume effect. Further study at multiple volumes and pion masses is 
necessary to disentangle the chiral and volume corrections. 

B. Charged pion and kaon 

We utilize the conventional effective mass plot in order to display the non-standard be- 
havior of charged particle correlation functions. In Fig. [31 we display effective mass plots 
for the charged pion and charged kaon. In non-vanishing fields, correlators exhibit a clear 
rise in the effective mass, Eq. (II ip. with respect to time. The need for a fully relativistic 
treatment of the two-point function is also evident from the figure as effective-mass shifts 
are on the order of the rest mass. 

Fits to the correlation functions of charged particles have been shown in the effective 
mass plots, Fig. [3j We fit the charged particle correlation functions using contributions 
from two states, as in Eq. dHJ) Although the amplitudes of the two states, Z{£) and Z'(£), 
enter the fit function linearly, we have not utilized variable projection due to the increased 
computational time needed to perform the fits. In zero field, we augmented the correlation 
function with backwards propagating contributions to the two states. In non-vanishing elec- 
tric fields, however, we found that backwards propagating charged particles make negligible 
contributions to the correlation functions. For a 1% effect due to a backwards propagating 
state, one must go beyond t = 40 for the n = 1 field strength, and to even larger times 
in stronger fields. Consequently we ignore backwards propagation in all but the zero-field 
case. Carrying out time-correlated fits on the bootstrap ensemble, we arrive at an ensemble 
of rest energies, {'£{(£)}, for the ground state. At this point, the analysis parallels that of 
the neutral particles. Fits to the energy function are carried out on the bootstrap ensemble 
using Eq. ffl~2]) producing the mass, polarizability, and quartic coupling. These extracted 
parameters are then averaged over the bootstrap ensemble. Their uncertainties arise from 
both fitting and bootstrapping, which we have added in quadrature. 

Extracted values of rest energies and fit parameters have been tabulated for the charged 
pion and kaon in Table ILT1 In performing these fits, we used the fit window 8 < t < 29. By 
comparing fits on adjacent time windows, we can estimate the systematic due to the choice in 



11 



n + : n=l 



K + : n=l 




C 0.20 
— 

§ 0.15 

o.io 





T£ 0.20 
"3 0.15 



i£ 0.20 
"3 0.15 



i£ 0.20 
<3 0.15 



15 20 
t/fl, 

# + : n=2 




5 10 15 20 25 30 

t/a, 
K + : n=3 





FIG. 3: Charged pion and kaon effective mass plots along with two-state fits to the correlation 
function using Eq. ([8]). The bands represent the uncertainty from the extracted ground state rest 
energy, which is also plotted separately as a flat band. Using the fit parameters Z{£) and E(£), 
we form the ground state correlation function which is depicted by the intermediate (gray) bands. 
Over long times, the effective mass should asymptotically approach these bands. We omit the zero 
field plots because they are quite similar to their uncharged counterparts. 
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0.10 




5 10 15 20 25 30 5 10 15 20 25 30 



t / a, t/a, 

FIG. 4: Effective energy plots for the charged pion and kaon two-point functions. Fits to the 
correlation functions are also shown. Values for n correspond to the strength of the quantized 
electric field £ , given in Eq. ([3]). We have omitted the n = plots, as they are just the ordinary 
effective mass plots. 

fit window. We find a large spread in the extracted value of charged particle polarizabilities, 
and consequently a comparatively large systematic uncertainty due to the fit window. Rest 
energies are particularly sensitive to the fit window as the field strength increases. 

While fits to charged particle correlation functions appear to describe the data well when 
displayed in terms of the effective mass, a further tool can be used to more clearly present 
these fits. This tool, moreover, aids in the determination of appropriate fit windows, and 
we refer to it as the effective energy plot. The effective energy, just like the effective mass, 
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FIG. 5: Plots of the electric field dependence of the extracted rest energies for charged mesons. 

is produced by considering the correlation function at successive times. The relativistic 
propagator for a charged particle in Eq. (jHJ) depends on the time, the electric field, and rest 
energy, D = D[t,E(£),S), albeit through a complicated one- dimensional integral. Given 
numerical data for the correlation function, g(t, S), we can successively solve 6 for the effective 
energy in time by considering the ratio 

D(t + l,E efi ,£) = git + 1,8) 
D(t,E e{i ,S) g(t,£) ' 1 ] 

with the value of the electric field, £, as input. This produces the effective energy as a 
function of time, E e g(t). Effective energy plots for the charged pion and kaon are shown in 
Fig. HI The effective energy should plateau over long times to the rest energy of the charged 
particle. From the figure, however, we see that contributions from the first excited state 
linger, and plateaus are not quite reached before the noise grows substantially. Nonetheless, 
we clearly see behavior reminiscent of the neutral particle effective mass plots in Fig. [TJ This 
confirms that Eq. properly describes the correlation function of a charged particle in an 
electric field. 

Finally in Fig. \5\ we plot the electric field dependence of the extracted rest energies of 
the charged pion and kaon. There is striking non-monotonic behavior which indicates the 
presence of quartic and perhaps higher-order terms in the field strength. We can make a 



6 Because the effective energy is deduced from the non- linear relation in Eq. (fT4|) . there is no guarantee 
a solution exists. Ensembles for which no solution can be found at a given time are dropped from the 
bootstrap. This only affected error bars the n — 4 effective energy plot for the ir + , and only for t > 24, 
where on average 5 bootstraps were dropped. 
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brief comparison with chiral perturbation theory. The size of the extracted polarizabilities 
is consistent with naive expectations, i.e. positive and on the order of 10 _4 f m 3 in physical 
units. 



V. CONCLUSION 



In this work, we have employed constant electric fields on a periodic lattice to investigate 
meson electric polarizabilities. Sizes of current-day lattices allow the utilization of properly 
quantized values of the electric field that lead to perturbative shifts in hadron energies. To 
test our setup, we have shown that the neutral pion (connected part) and kaon polarizabilities 
can be extracted from lattice QCD by measuring their energies as a function of the applied 
electric field strength. Furthermore, we have investigated the charged pion and charged kaon 
polarizabilities, for which simple spectroscopy is of no avail. Using the relativistic charged 
particle propagator in the presence of an electric field, we fit lattice two-point functions and 
extract rest energies of charged pions and kaons. Using effective energy plots, we showed 
that, despite non-standard behavior for the correlation function, rest energies of charged 
particles show behavior similar to the effective masses of neutral particles in electric fields. 
Charged meson polarizabilities were then extracted from the behavior of the rest energy as 
a function of the electric field. Resulting electric polarizabilities have comparatively large 
uncertainties due predominantly to two sources. With our current analysis method, the 
choice of time window gives a larger than expected systematic uncertainty. Global fits, 
that are correlated in both time and electric field strength, can address such systematic 
error. Secondly higher-order terms in the weak field expansion of charged particle rest 
energies appear to be very important, prompting future study on larger lattices on which 
the quantized field strengths are smaller. We hope that further refinements to the fitting 
procedure, additional data at different volumes and pion masses will remove the largest 
systematic effects, and ultimately bring lattice QCD in contact with experimental data for 
polarizabilities. Ultimately we will also use sea quarks that couple to the background fields. 
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APPENDIX: NON-UNIFORM FIELDS 



Although we employ quantized field strengths, Eq. ([3]), with a proper treatment of the 
boundary flux, Eq. (j3J), we have additionally explored the effect of non-quantized fields on 
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FIG. 6: Comparison of background field implementation. On the left, we fix the field strength and 
change the source location using Eq. ([2]) to implement the background field. On the right, we fix 
the source time and change the field strength while using Eq. (Tjrj). 



particle correlators. For this study, we use isotropic 24 3 x 64 lattices, the details of which 



are presented in [21|. We summarize our findings here. 

First we consider the naive implementation of the external field using Eq. (j2J) and the field 
value corresponding to n = 3 in Eq. ([3]). In the continuum limit, the spike in the boundary 
flux contracts to a point and the field becomes uniform. To test the uniformity of the field 
at finite lattice spacing, we look at the (connected) neutral pion two-point function. If we 
take the source time at t STC = 0, then the effective mass exhibits a plateau around t = 12, as 
shown in Fig. O On the other hand, if we take the source time at t src = 52, then the plateau 
would set in as the pion wraps around the time boundary. The correlation function shows 
striking evidence for the spike in the electric field from boundary flux. Notice in plotting we 
have translated the latter correlation function forward by 12 units in time. 

Next we consider the proper implementation of the external field on a torus using Eq. (j3j). 
We again consider the (connected) neutral pion two-point function. Fixing the source time 
we plot in Fig. E]the resulting effective mass for two values of the field strength, 
n = 3 and n = e = 2.71828 ... . We translate both correlation functions forward by 12 
units in time. For n = e, the effect of boundary flux has been mitigated (roughly by a factor 
of ten), but, leads to easily measurable shifts in the particle energy. The quantized value, 
n — 3, exhibits a plateau as the field is uniform across the time boundary. 
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